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Abstract 

We consider the effect of surface roughness on soHd-soUd contact in a Stokes flow. Various 
models for the roughness are considered, and a unified methodology is given to derive the corre- 
sponding asymptotics of the drag force. In this way, we recover and clarify the various expressions 
that can be found in the litterature. 



1 Introduction 

The dynamics of solid particles in a viscous fluid is crucial to many phenomena, such as blood flow, 
sedimentation or filtration. The drag force exerted by the fluid on the solids plays of course a central 
role in this dynamics. It has been the matter of many studies. The first ones focused on the dynamics 
of a rigid sphere near a plane wall, that moves in a Stokes flow under no-slip conditions: we refer 
to the pioneering works |[5] |2T] |7] |22]. The main conclusion of these works is that the drag force 
is inversely proportional to the distance h = h{t) between the sphere and the plane at time t. The 
reduced ordinary differential equation that governs the movement of the sphere is then of the type: 
h + h/h = f, which prevents collision between the sphere and the wall in finite time. We quote that 
this striking conclusion holds for any value of the fluid viscosity and of the sphere density. Moreover, 
it is still valid for arbitrary solids with smooth surfaces, and it is still valid within an unsteady Navier- 
Stokes flow (see 11121 ). 

This theoretical no-collision result, that goes against Archimedes' principle, is clearly unrealistic 
at the scale of macroscopic solids. Even at microscopic scales, "dry collisions" have been clearly 
recognized. Therefore, many articles have tried to identify the flaw of the previous modelling, in 
order to circumvent the paradox. Among possible flaws that have been suggested one can mention: 

• The rigidity assumption. Elasticity, even weak could allow for solid contact: see ||8l. 

• The no-slip condition, that is no longer valid when the distance between the solids is of the 
order of the mean free path of the fluid particles: see llT4l . 

• the incompressibility assumption : see [2]. 
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We shall focus here on another very popular explanation for the no-collision paradox: roughness. The 
basic idea is that nothing is as smooth as a plane or a sphere: iixegularity of the surface can thus affect 
the fluid-solid interaction. This credo has led to many experimental and theoretical studies, focusing 
on roughness-induced effects on drag forces (||23] S [19]). Such studies will be discussed in the core 
of the paper. 

We quote that the interest in roughness issues has been renewed these last years, notably in con- 
nection to microfluidics. Indeed, it has been recognized that the classical no-slip boundary condition, 
which is relevant at the macroscopic scale, may fail at the micro- or nanoscale. This happens for in- 
stance for some corrugated hydrophobic surfaces, which trap gas bubbles in their humps and generate 
in this way some substantial slip. More generally, to determine the appropriate boundary condition at 
a rough surface is a matter of current debate. In this context, if one has theoretical formulas that ex- 
press how the drag force depends on the "rough" boundary conditions, one may check experimentally 
through the force measurement what the right boundary condition is. This interesting point of view is 
for instance developed in |[r7ll24l . 

The aim of this paper is to investigate mathematically and in a unified way the relation between 
the roughness and the drag force. Namely, we study the evolution with time t of a rough solid S{t), 
falling towards a rough wall P in a Stokes flow. We assume for simplicity that the solid moves by 
translating along the vertical axis r = 0, where (r, 0, z) are cylindrical coordinates. We shall comment 
on this simplification later on. Various models for the roughness are to be considered. In all models, 
the moving solid is described at time t by S{t) = h{t) + 5 for a fixed S. We assume that S has its 
lower tip at r = 0, and that in the vicinity of its lower tip, its surface is described by: 

z = ls{r), r<ro, 6 £ {0,2tt) 

for some < 1 and some Lipschitz function 75 with 75(0) = 0, 75 > 0. Notice that the solid 
velocity is given by h{t) Cz- Similarly, the wall P is described in cartesian coordinates (x, y, z), by 

z = lp{x,y), (2;,y)G]R^, 

for some Lipschitz function 7p, with 7p(0, 0) = 0, 7p < 0. Accordingly, we denote the fluid domain 

F{t) := {yi= {x,y,z), x^S{t), z>'yp{x,y)}. 

lfu = u{t,x) = {ux{t,x),Uy(t, x),Uz(t,x)) andp = p(t,x) stand for the fluid velocity and pressure, 
the steady Stokes equations read 

-Aii + Vp = 0, divn = 0, t>0, xGF(t). (1.1) 

We neglect gravity, as it plays no role in the discussion. Our goal is to study the force on the sphere, 
that is 

Fd{t) := / {2D{u)n-pn)da ■ ez- (1.2) 
JdS(t) 

The notations n and D{u) refer to the nomial vector pointing outside the fluid domain and the sym- 
metric part of the gradient respectively. 

In order to determine Fd{t), one needs to specify the boundary conditions at the solid surface and 
at the plane. In all our models for roughness, such conditions have the following general form: 

[u - h{t) ez) ■ n\Qs(^t) =0, {u - h{t) e^) x n\QS(^t) = -2/3s [D{u)n] x n\QS(^t) (1-3) 
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and 

u-n\p = 0, uxn\p = —2l3p[D{u)n]xn\p (1.4) 

where Ps^Pp S [0,+oo). These ai^e boundary conditions of Navier type, the constants /3p and /3s 
being the slip lengths. Of special importance is the case = I^P = 0, which corresponds to the 
no-slip condition. 

We model the roughness in three different ways: 

1. through a lack of differentiability. Namely, we consider a solid S which is axisymmetric around 
{r = 0}, and satisfies 

75.(r) := 1 - - r2 +er^+", a G [0, 1), r < tq. 

This means that the solid surface is locally a smooth sphere z = I — \/l — r"^, perturbed by a less 
regular "rough profile" of amplitude e. For a = 0, this profile is a spike, which has Lipschitz 
regularity. For a > 0, the profile is differentiable, with a Holder derivative. For simplicity, we 
do not consider any roughness on the wall, and take the classical no-slip boundary conditions: 

7P = 0,^p = (3s = 0. 

2. through a slip condition. We consider the case of a ball S, of radius 1, falling vertically above 
a plane wall, with positive slip coefficients: 

7s(r) = 1 - \/l -r2, 7p = 0, (3s, f3p > 0. 

Let us stress that such modelling of the roughness by the addition of (small) slip is commonly 
used. It is well-accepted in the context of rough hydrophobic surfaces lH, and a topic of debate 
in the context of hydrophilic ones cf 10711241 . 

3. through a small parameter. Namely, the roughness is modelled through a small amplitude, high 
frequency perturbation of a plane wall. That means P is described by the equation 

z = lp{x,y) := e7(x/e,y/e), e <C 1 

for some periodic and smooth non-positive function ^{X,Y), with 7(0,0) = 0. In parallel, 
we assume no roughness on the solid surface (75 (r) = 1 — \/l — r^), and the classical no-slip 
conditions: /3s = /3p = 0. 

Note that if we take the parameters e, /3s and /3p to be zero in the previous models, we are back 
to the classical situation of a curved and smooth solid falling towards a plane wall. The whole point 
is to derive the next order terms that are involved in the expression of Td- Note also that, in view of 
our models, the assumption that the solid translates along r = is natural. For the first two models, 
the whole geometry is axisymmetric. For the third one, one can consider rough walls P"^ that ai^e 
symmetric with respect to x and y. In all these configurations, if the initial velocity field of the solid 
is along r = 0, both the geometry and the Stokes flow inherit strong symmetry properties, forcing the 
velocity field of the solid to be along r = for all time. 

The ambition of this paper is to provide a rigorous and general methodology to derive the drag 
term Td, in the regime of small distance h between the solid and the wall. This methodology, which 
relies on the calculus of variations, will be explained in section |2] Then, in section |3j it will be applied 
to our first two models of rough surfaces. In this way, we will extend results from former formal 
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computations, notably those in |[T4l[T9l . In the last section|4j we will turn to the third model of a small 
amplitude and high frequency boundaiy. This model is of particular interest, as it is connected to the 
phenomenon of apparent slip, which is a topic of current interest in fluid mechanics, see [18]. We 
will notably discuss the introduction of an effective sUp length as a modelling for hydrophilic rough 
surfaces. 

2 Methodology for drag derivation 

We present in this section a general approach to the derivation of the drag force Fd{t) on the solid 
sphere S{t). We first remark that the geometric configuration at time t is entirely characterized by 
the distance h{t) between the lower tip of the solid and the origin x = 0. Thus, we can rewrite 

S{t) = Sh(t),F{t) = Fh(^t), with the family {Sh, Fh)h satisfying 

Sh = h + S, F/i = |x, x^Sh, z>7p(x,y)|. 

Moreover, considering the linear- Stokes equation (II. Il l and boundary conditions (ll.3l )- (ll.41 i. we can 
write u{t, x) = h{t)uh(^t){x) and p{t, x) = h{t)pm^{x) where Uh,Ph satisfy the steady problem 

- Auh + Vph = 0, divu,, = 0, xeFh (2.1) 

together with the boundary conditions 

{uh - ez) ■ n\ds^ = 0, {uh - e^) x n\Qs^ = -2f3s[D{uh)n] x n\Qs^ (2.2) 

and 

Uh ■ fAp = 0, u/i X n\p = — 2/3p [D{uh)n] x njgp (2.3) 
Accordingly, we can write 

Tdit) = h{t)J^h{t), ■■= / {2D{uh)n-phn)da ■ Cz- 

JdSh 

The problem is to determine the behaviour of Fh in the limit /i — )• 0. Our method to address this 
problem has three main steps: 

1. In a first step, we express the drag Fh as the minimum of some energy functional. One can do 
it using the variational interpretation of (12. U - (l2.21 i- (l2.3l ). It allows to identify for all our models 
of roughness an energy functional 8h and a set of "admissible fields" Ah such that 

Fh = min £h{u). 
u&Ah 

The expUcit definitions of £h and Ah will be given at the end of this section. 

2. In a second step, we rely on the minimization problem introduced in Step 1 to find an accurate 
lower bound for Fh- Namely, we choose some appropriate energy functional £h < £h and 
some appropriate set of admissible fields Ah D Ah for which we can compute explicitly the 
minimimum and corresponding minimizer u. In this way, we get 

£h{u) = min £h{u) < min £hiu) = Fh 

which yields a lower bound. Of course, the relaxed functional £h and admissible set Ah must 
remain close enough to the original ones, in order for this lower bound to be accurate. We will 
make them explicit for our various roughness models later on. 
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3. In a third step, we choose some appropriate field u £ Ah so that 



= min Shiu) < £h{u) 

provides an accurate upper bound for the drag (that is with the same type of behaviour as the 
lower one). In many cases, as will be seen later on, the minimizer u G Ah of the second step 
generally belongs to the original set of admissible fields Ah, or at least can be slightly modified 
to belong to Ah- Thus, one can take in general u ^ u. 

Our goal in the present paper is to apply this methodology to have a better understanding of 
roughness effects. In this section, we cariy out step 1, that is the formulation of the drag in temis 
of some minimization problem. This step is very general, and independent of the roughness issues. 
In the next sections, when turning to step 2 and step 3, each roughness model will of course require 
specific calculations. 

To link the drag to an extremum problem, we must distinguish between the case of no slip (/3s = 
/9p = 0) and the case of non-zero slip {/3s > 0, /3p > 0). 

• In the case of no-slip, the divergence free-condition implies 

\Vuh\^ = 2 [ \D{uh)\^. 



Hence, multiplying the Stokes equation (12.1b by Uh and integrating over the fluid domain Fh, 
we obtain by Stokes formula 

2/ \D{uh)\'^ = / {2D{uh)n- phn)-Uhda = / {2D{uh)n - phu) da ■ = Fh- 
JFh JdShUP JdSh 

Moreover, we know that equation (12. Il l (together with the boundary conditions (I2.21 i- (I2.3I )) is 
the Euler equation of a minimization problem. Namely, 

j |Vuh|^ = min|y" |Vup, u£Hl^^{Fh), V • u = 0, u\p = u\s^ = 

(We remind that the Sobolev space H^^^ is the space of fields u that are locally square integrable, 
with distributional derivative Vn also locally square integrable). Q 

Indeed, if u has the properties mentioned above, then u—Uhi^ zero along the boundary dSh^P- 
So, multiplying (11.11 ) hy u — Uh and integrating by paits, we end up with 

/ \Vuh\^ < I Vuh-.Vu < ( [ \Vuh\A ' ( [ \VuA ' , 

using the Cauchy-Schwarz inequality. The characterization of Uh follows, and eventually yields 
thatJX = min„g^^ with 

Ehiu) := |Vu|2, Ah ■-= |n G Hl^{Fh), V-u = 0, u\p = 0, = e^} (2-4) 



' As no abstract theory is needed in ttie remainder of the article, such mathematical details can be skipped without harm. 
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• In the case of positive slip lengths /J^, /3p, the computation is slightly different: multiplying the 
Stokes equation hy — u, u E Ah, we obtain after integrating by parts: 



/ 2D{uh) : D{uh-u)+ ^ [ (K 
Jf. PS JdSh 



Cz) X n) • {{uh — u) X n) da 



+ — / X n) ■ {{uh. — u) X n) da = 0. 
IP 



/5p 



In order to recover full gradients instead of symmetric gradients, we proceed as follows. On one 
hand, by standard identities of differential geometry (see for instance ||6] Lemma I, p. 233]), 
we have ^ ^ 

D{v)n X n = -^jdnV x n + -^v x n at dSh 

for any smooth v satisfying v • n = at dSh. The last term at the r.h.s is connected to the 
curvature of dS, which is simply 1 by our choice of S. Similarly, 

D{y)n X n = ^'^nV x n at P 

for any smooth v satisfying v ■ n = at P. On the other hand, writing 

Av = div (Vv), respectively Av = 2div {D{v)), 

and integrating by parts, we get that for any smooth v in Ah and any smooth w satisfying 
-u; • n = at dSh U P: 



Av ■w = —l Vv : Vw + / {dn{v — Cz) x n) ■ {w x n) + / x n) ■ {w x n), 

Fh JFh JdSh Jp 

respectively 

Av-w = — / 2D{v) : D{w)+ {2D{v — ez)n x n)-{wxn)+ I {2D{v)n x n)-{wxn). 
Fh JPh JdSh Jp 

Combining the previous identities, we get 



2D{v) : D{w) = / Vv:Vw+ / {{v - Cz) x n) ■ {w x n). 
Fh JFh JdSh 

We take v = Uh, w = u — Uh, and inject this last equality into (12.51) to obtain 

/ |Vn/j|2 + ( J- + 1 ) / \[uh - ez) X n\'^ da + ^ I \uh X n\'^ da = 1 Vuh'.Vu 
JPh \PS / JdSh PP Jp JPh 

+ ( ^ + 1 ) / {{uh - ez) X n) ■ {{u - ez) X n) da + ^ I {uh x n) ■ {u x n)da. 
\ps J JdSh pp Jp 

Use of the Cauchy-Schwartz inequality and of the Young inequality \fab < ^(a + 6) leaves us 
with 

/ \Vuh\^ + {^ + 1] I \{uh - ez) X da + ^ I \uh X da 
JFh \PS J JdSh PP Jp 

< I iVii]^ + ( + 1 ) / - ez) X n\^ da + I \ux n\'^ da. 

JPh \PS J JdSh PP Jp 
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Thus, we have this time that Jtj = inin„g^,^ with 

£hiu) := I |Viip+ (-^ + 1) / \{u-e^) xn\^ + ^ I \uxn\^ 



Ps J JaSh Pp Jp 



Ah := |u G Hl^{Fh), V • ti = 0, u-n\p = 0, (u-Cz)- n|as,^ = o|. 



(2.6) 



We note that contrary to the no-shp case, only the impermeabiUty condition is included in the def- 
inition of the space Ah- It can be shown that the Euler equation for the latter minimizing problem 
includes the boundary conditions (I2.2I )- (I2.3I ) on the tangential part of the velocity-field by standard 
integration by parts as in the no-slip case. For brevity, we shall replace the coefficient 1/ f3s + I by 
1/ (3s in what follows. This means that we shall include curvature effects in the slip coefficient. The 
characterization of the drag through energy functionals (12.41) and (12.61) will be applied to our first two 
roughness models in the next section. 



3 Application to various roughness models 

In this section, we detail the steps 2 and 3 of our methodology, both in the case of a non-smooth 
boundary (model 1) and in the case of slip boundary conditions (model 2). 



3.1 The case of non-smooth soUds 

As emphasized in the introduction, we consider here the case of an axisymmetric solid S, whose 
boundary is described neai^ its lower tip by 

75(r) = 1 - v^l - r2 +er^+", aG[0, 1], r < tq. 

The wall is flat, and no slip conditions are imposed at all boundaries. The drag is given by 

Th = min £hiu), 

u&Ah 

with the energy £h and the set of admissible fields Ah given in (12.41 ). 

As the fluid domain Fh is invariant by rotations around e^, much can be said about the minimizer 
u = Uh- Indeed, for any rotation Rq around e^, RgUhR-e still belongs to Ah, and has the same 
energy as Uh- Uniqueness of this minimizer yields 

R0Uh{R-ex) = Uh{x), VxgF/i. (3.1) 

This means that Uh has the following structure: 

Uh = Uh,rir,z)er + Uh,9{r,z)ee + Uh,zir,z)ez, 

where (r, 0, z), resp. (e^, eg, e^) are the cylindrical coordinates, resp. the cylindrical vector basis. One 
then remarks that Vh = Uh^ri^, z)er + Uh^zi'r, z)ez still belongs to Ah, with £h{vh) < £h{uh)- Again, 
by uniqueness of the minimizer, we get Uh = Vh and Uhfi = 0. Thus, the divergence free condition 
resumes to ^ 

-dr{rUh^r) + dzUh,z = 0. 
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Together with the boundary condition Uh^.z{r, 0) = 0, it leads to 

Uh = -dzct)er + -dr{r4>)ez, (3.2) 
r 

with sti"eamfunction (f){r, z) := — Uh^ri''", z') dz'. The boundary conditions on are 

d,<l)\as = d,<l)\p = 0, dr{r^)\as = r, </.|p = 0. (3.3) 

Thus, we can without restriction include these last conditions in the set of admissible fields: instead 
of the original definition in (I2.6l l. we take 

Ah ■= |ii G Hlg^{Fh), u = —dz4>er + -(9r for some cj) satisfying (I3.3I )|. 

We quote that the boundary conditions on at dS yield 

5^(/)(r, /i + 75(r)) = 0, dr{r(t)){r,h + -is{r)) = r, r < r^. 

They imply in turn that (/)(r, h + "fs{r)) = | + ^ for some constant c. As (/)(r, z) = — Jq Ur{r, z') dz' 
is regular enough neai^ r = 0, we deduce c = 0. Eventually 

r 

dz(l){r,h + -fs{r)) =0, H^, h + jsir)) = -, r < ro, 

^ (3.4) 

dz(j){r, 0) = 0, (/>(r, 0) = 0, r < tq. 

From there, we obtain an accurate lower bound as follows. Noticing that 

|Vup = \drz<P^\ + \dzHr\^ + \dzz(l>\^ + \dr[dr{rcl))/rf + |a„(r0)/r|2. 

we anticipate that in the limit of small h, most of the energy will come from a neighborhood of the 
lower tip of the sphere := {r < r^, < z < h + 75(r)}, and will be due to the z derivatives 
of the stream function (p. Accordingly, we introduce the following relaxed minimizing set and energy 
functional: 



Ah ■= ^uGH^{F^), u = —dz(t>er + -dr{r((>)ez for some 6 satisfying (13.41) |. 



Sh := / \dzUr\^ = / 



(3.5) 



From the Euler equation 5f = and the boundary conditions (13.41 ). it follows easily that the 
latter minimum is realized with 

<^/,(r,z) = ^$(— ^— V where m=t\'i-2t), ViG[0,l] 
and has for value: 

Th = 6TT J— ^-"3 . (3.6) 

Jo ih + ls{r)r 
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We emphasize that this formula is general for no-slip boundary conditions. It does not require any 
special assumption on the solid surface. In the case 7s'(r) = 1 — \/l — r'^ + er^+", our lower bound 
satisfies 



r^dr 



(/i + ^ + eri+° + 0(r4))3 



^X^e/i"^) + 0{J{eh^,h)) + 0{l) 



where 



X(/3) := 



s^ds 



s ds 



(3.7) 



— , J{p,h) : = 

/O (l + ^+^sl+")3 Jo (1 + f + /3sl+-)4 

The computation of the asymptotic behaviours of I and J' is detailed in Appendix [A] It yields the 
following results: 



When /3 ^ 1, we obtain : 



1 



with an explicit constant Aq, given in the appendix. 
• When /3 » 1, we have: 



+ 0(/32), J{/3,h) = 0(|ln(/i)|) 



(3.8) 



' / 1 



9 ln(/3) 



/33 



for a > 1/3, 
for a = 1/3, 
for a < 1/3, 



(3.9) 



where the value of is also provided in the appendix. As regards the remainder, we have the 
following bound: 

J(/3,/i)=0(|ln(/3) + i^ln(/i)|) 



Back to the drag force, (13.81 ) and (|3.91 l yield the following lower bound: for /3 = e/i 2 ^1 



h + XaSh 2 



— (l + 0(/3)) + 0(|ln(/i)|) 



(3.10) 



g-l 

and for (3 = eh 2 > 1 



l~3a 



^ (I + /3— ) +0(|ln(/i)|), for a > 1/3, 



97r|ln(/i)| 

2 £3 



+ 



ln(e) 



6^/^° /l + /3^)+0(|ln(e)|), 



2 



for a = 1/3, 



for a < 1/3, 



(3.11) 
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Note that the expression given in the case a < 1/3 only matters when e <^ 1 (otherwise, one can just 
retain that = 0(1)). 

This concludes our study of a lower bound for the drag. Such bound is accurate, as we can with 
minor modifications obtain a similar upper bound. Indeed, it is an easy exercise to find a regular 
stream function (ph = (l)h{r, z) defined on Fh, equalling say on F^r\{r < ro/2}, such that 

% = V X {(j)heg) G Ah, 

and such that 



/ \dUh\^ = [ \dl~4>h\^ + 0(1 



uniformly in h and e. We quote that the remainder term is uniformly bounded, because no singularity 
is created outside of the contact zone (that is outside a vicinity of r = 0). We refer the reader to 
lfT2lfTn[T3l for more details in the case e = 0. Hence, we have: 

J^h < £h{uh) 

= [ [\drz^h\^ + \d,^h/r\^ + \dr[dr{r^h)/r]f + \drz{r^h)/r\^]+0{l) 

JFh ^ 

'|a.,0,.|2 + |a,0fc/r|2 + |c>,[a,(r0^)/r]|2 + |a,,(r0;,)/r|2l + 0(1) 



Th + 



The computation of the integral terms at the r.h.s. follows the lines of Appendix |A] It yields some 
0(min(| ln(/i)|, | ln(e)|)) error term. The main reason for these integrals to be lower order terms is 
that in the "curved" contact zone, the typical lengthscales in z and r are respectively h and \/7i, so 
that z-derivatives ai^e more singular than r-derivatives. Eventually, 

J^h = h + 0(min(|ln(e)|,|ln(/i)|)) 

We stress that this modelling of the roughness solves the famous no-collision paradox discussed in the 
introduction. Indeed, as h goes to zero for a given e, (3 goes to infinity and the roughness effect yields 
a drag force Th which is always bounded by for some 7 < 1. In particular, it is weaker than in 

the smooth case, and the solid dynamics, which is governed by the o.d.e. 

h + hTh = 

allows for h to cancel in finite time. 



3.2 The case of slip boundary conditions 

We turn in this pai^agraph to our second model, in which roughness is involved through slip coeffi- 
cients. We want to have a close approximation of J7, = min_4,^ £h, where this time Ah and £h are 
defined in (12.61) . We still have a rotational invariance in this case, so that we can again reduce Ah by 
restricting to velocity fields of the type 

u = -dz4>er + -dr[rcl)\ez. 
r 

The impermeability condition at P yields again 

<^(r,0) = 0, VrG(0,ro). (3.12) 
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As regards dS, we have for all r £ (0, tq) 
1 



e. ■ n 



u ■ n 



1 

Vl + l7^(r)P 



so that the impermeability condition leads to 



lsir)dz(j){r, h + 75(r)) + -dr{rcl)){s, h + -fs{r)) 

r 

--^[r(/>(r,/i + 75(r))], 



(3.13) 



(r, /i + 7s(r)) = -, VrG(0,ro). 



(3.14) 



Accordingly, we inti'oduce the relaxed set 



Ah := |n G H^{F^), u = -d;,(j)er + -dr[r(l)] e^, <p satisfying (I3l2]) and (ITT4]) | . 

with defined in the previous section. We then need to define the relaxed energy S^. As in the 
previous section, we shall keep only dzz4' the gradient terms. But we shall not change the boundary 
integrals involved in (12.61 ). Therefore, we compute: 



on P, u X n = dz 



• ondSn{r < ro}, because of (l37T3l) - (l3l4l) . (u - e^) x n = y^l + \'y'g{r)\^dz4> eg. 
Hence, we introduce the approximate energy 



F°n{r<ro} 



|c^zz0('^> z)\'^rdrd9dz 



+ 271 



(3s 



|9,</>(r,/i + 7s(r))|2 + — |9,0(r,O)|2 
PP 



rdr. 



The corresponding mimimization problem is easy, because it amounts to find, for each value of r < ro, 
the minimizer of the functional 



£h{r) 



/i+7s(r) 



hl{z)\'^dz + 



1 



\4>'^{h + ^s{r))? + ^m^)? 
PS PP 



over functions c^j. = (f)r {z) satisfying the inhomogeneous Dirichlet conditions 



(/>,(0)=0, Mh + ls{r)) = -. 

This is a one-dimensional minimization problem, with Euler equation 
Dirichlet conditions, plus Robin type condition on 4>[:. 



0, endowed with above 



^':.{h + 75(r)) + + 75(r)) = 0, 



(is 



<^';(o) - -^0;(o) = 0. 

PP 
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After a few computations, the minimum of 8^ is obtained for 

r ^ f z 



Oh 



(r, z) = r 



2 \' h + 75(r) 

where $(r, t) is the polynomial of degree 3 in t given by 

^/ ^^ _ 2 {as + asap + ap) ^3 _^ 3 (2 + gg) gp ^2 

' 12 + 4 (05 + ap) + Qp as 12 + 4 (a^ + ap) + ap a^ 

6 (2 + as) 

'^1 



(3.15) 



12 + 4 (as + ap) + ap as 
where: 

{^ + \is{r)?fHh + ls{r)) {h + js{r)) 
as = as{r) := , ap = ap[r) := . 

PS PP 

Note that the coefficients of $ are uniformly bounded in as, ap, that is in r < ro, /3p, 13s, h. In the 
limiting case /3s = /3p = (no-slip limit), we obtain formally ^{t) = —2t^ + 3i^, in agreement with 
the computations of the previous section. 

We now turn to the lower bound 



J^h = mill Sh 

Ah 



fro r /"l 

27r/ / \dtt^{r,s)fds + as\dt^{r,l)\^ +ap\dt^{ 
Jo Uo 



r dr 



(/i + 7s(r))3 



We make the last integral more explicit by replacing by its value. We obtain 



= \ f\h{r) + h{r)) 



(/i + 7s(r))3 

where the integrands Ii and I2 ai^e given by 

12 (a| ap + 5 (a| ap + ap as) + 4 (a| + ap) + 20 as ap) 

Ii := 2 

(12 + 4 (as + ap) + as^p) 

144 (as + ap) 



(12 + 4 (as + ctp) + asap)'^ 



Note that /i and I2 are uniformly bounded in as, ap, that is in r < tq, /3p, /3s, ^- Thus, expanding 
7s, we obtain: 

= U^\h{r) + hir)) + O(J(0,/i)) 

ro ^ t/^ (3.16) 

(/i(r)+/2(r)) 1, + 0(|ln(/.)|) 



2^0 (/i + T)' 

where ^(0, h) was introduced in (I3.71 i and shown to be 0(| ln(/i) | ) . We must now distinguish between 
two cases, depending on the behaviour of h/ f3s and h/ jip: 
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1. Either h/ (5$ or /i//3p is of order 1 or lai^ger. Then, either as or ap is of order 1 or lai^ger. It 
follows that 

c < h{r) + hir) < C, 

for all r < vq, where the constants c, C are uniform with respect to all parameters. We then 
deduce from (13.161 ) that 

c' ~ C 

Note that in the limiting case /J^ = /3p = (no-slip limit), we obtain formally: 

Gvr 

as = ap = +00, h = 12, h = 0, Fh = — 

n 

recovering the classical result. We also emphasize that the regime considered here includes the 
case where one of the slip coefficients is zero. In paiticulai", the drag force is stronger than c' /h 
in such a case, preventing any collision. 

2. Both h/Ps and /i//3p are small. This case requires more care. We first notice that 

«P = ;-^(^ + y + 0(^')), as = j^(h + ^-^+0{r\h + r'))y (3.17) 
From there, for ro and h small enough, we get 

cJi(r) < /i(r) < CJi(r), J,{r) := i + 



where c, C > and 



1 + ap 1 + as 



1 /, r'^\ , . 1 /. 



Then, with the change of variable r = , we write 

Mr) 



+ 



1 V/3p (1 + ^^72) ^ V/35 (1 + \ ^ u^du 



hJo \^l + V/3p (1 + ^2/2) l + h//3s{l + uy2)J (l + tx2/2)3- 

In the regime of small /i//3p and h/f3s, we get that this last integral is o(l//3p + 1//35). Finally, 
all of this leads to 



- h{r) ^ = o(l//3p + I//35). 



(3.19) 



It now remains to evaluate the contribution of I2, which will yield the leading behaviour of Fh- 
The use of (13.17b gives first 



vr 



^ ^ r-^dr vr f'^ 144 (05 + ap) r^dr 

3 ~ 



{h+'iY 2 Jo {12 + A {as + ap) + asapY {h + '-^y 

+ 0(l//3p + I//35). 
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Then, straightforward manipulations show that 



{1 + ciasY {1 + ciapYJ (12 + 4 (05 + ap) + agap)^ 

V(l + C2as)2 (1 + C2ap)^y 

for some ci,C2 > 0. As seen in the treatment of Ii, the 0(Ji(r)) term will only contribute 
to the drag through a o(l//?p + 1/Ps) term. The main contribution of I2 to the drag will be 
governed by 

TT f as ^ ap \ r^dr 



2 Jo \il + cas)^ (l + cap)V {h + 

TT f^. ( h/ps (1 + nV2) ^ /i//3p (1 + ^2/2) ^ -u^dn 



2/1 7o \{l + ch/Ps{l + uy2))^ (1 + cV/3p (1 + ^^72) )V (1 + ^^)3 
""l^ii {l + ch/Psx)^x^ ^ J^Ji {l + ch/Ppxfx^j 



vr + I ln(/i)| + 0(l//3p + l//?^) 



through standard manipulations. It yields eventually 



^ rh{r) 77^^ = ^ + ^) 1 1^(^)1 + 0{1/Pp + (3.20) 
2 Jo (/i+V) VPs ppj 



Combining (13.161 ). (13.191 ) and (13. 201 ). we end up with the following lower bound for the drag: 

-^/^ = vr (-^ + I ln(/i)| + 0(l//3p + l/Ps) + 0(1 ln(/i)|). 

This lower bound is similar to the one derived by L.M. Hocking (see HH). 

This concludes our study of a lower bound for the drag. Hence, it remains to obtain a similar 
upper bound. One could develop the same approach as in the previous section. Namely, one could 
look for some suitable extension c^/^ of cj)}^, with similar behaviour for its energy. However, due to 
the elaborate expression (13.151 ). this would lead to tedious computations. We overcome this technical 
difficulty as follows: 

1. When h/ Pp or h/ Ps is of order 1 or larger, we take u/i = V x {<j)hee), with the "no-slip" 
streamf unction cj)}^ built in the previous section. We obtain with this choice some 0{l/h) upper 
bound as expected. 
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2. When /i//3p and h/ (5$ ai^e small, a good way to recover the right asymptotic behaviour is to set 

% := V X [(t)h{r,z)ee] with 

4>h{r,z) := ( r, 



2 V h + js{r) 
in F^, where 



1 1 \ t / ap as \ t"^ 



1 + ap 1 + 05/2 \l + ap 1 + as 

We extend then to the whole of Fh with a stream function having bounded gradients. Calcu- 
lations similar to the previous ones yield: 

£h{uh) = 7r(^^ + ^) |ln(/i)|+0(l//3p + l//3s) + Oi\lnh\) 

where we insist that the 0(| In is uniform with respect to /3p and Ps- In particular, in the 
realistic regime of small slip lenghts, we obtain the exact same leading behaviour for the lower 
and upper bounds. For the sake of brevity, we leave the details to the reader. 



4 The case of a corrugated wall 

In this section, we focus on the third model of roughness described in the introduction, in which the 
wall has a small amplitude and high frequency oscillation: namely, 

— , - j , e > 0, 7 = 7(X, Y) 1-periodic, 7 < 0, max 7 = 7(0, 0) = 0. 

We remind that the solid is assumed to be smooth, and that no-shp conditions hold both at the solid 
surface and the wall. We shall pay special attention to the regime e ^ /i <C 1, that is when the 
distance between the solid and the wall is much greater than the size of the roughness. 

Such roughness model with a small parameter is very popular, as it allows for multiscale analy- 
sis. This analysis has been notably performed in the context of wall laws. In this context, the idea 
is to replace the rough boundary by a flat one, and to impose there some good homogenized bound- 
ary condition, that expresses the mean effect of roughness. This homogenization problem has been 
considered by physicists since the early 90's, through numerics and explicit calculations for special 
geometries: see for instance Il20l . It has been adressed later on in some mathematical works, based 
on homogenization theory. We refer to |[T] [15] for periodic patterns of roughness, and to E [TO] for 
random roughness. The conclusion of these works is that, for small enough e, one can replace the 
oscillating boundary by the flat one {z = 0}, and impose there some Navier-type boundary condition: 

Uz = 0, {ux.Uy) = eBdz{ux,Uy) 

for some two by two positive matrix B, which is sometimes called the "mobility tensor". There has 
been a recent interest on qualitative properties of this tensor, for instance for shape optimization in 
microfluidics, cf |[T6ll . 

Another frequent idea is that a slip condition amounts to a no-slip condition at a shifted wall. 
Combining this idea with the previous one, some recent articles have suggested a drag force of the 
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type Th ~ h+pe some positive /?: see ifTTlfTQl . We will discuss this result in a rigorous manner 
here. 

First, one can use the methodology of section |2]to derive some lower and upper bounds. 

• As regai^ds the lower bound, let us show that 

Gtt 

J'h > 1 ^ + 0(|ln(/i + Ae)|), for A := -min7>0. 

n + Xe 



Indeed, we have 



Fh = min £h{u) 



where £hiu) and Ah are given by (12.41 ). Let us now define 

P^:={z = -eA}, := {x, x S/,, z > -eX}. 

Any field u of Ah can be extended by zero below the rough wall so that it can be seen as an 
element of the larger set 

Ah ■■= e hI^{F^), V-u = 0, u\px =0, u\os^ = ez'^. 

Then, obviously, 

Fh > mill £h{u). 

ueAh 

But the r.h.s of this inequality is exactly the drag force associated to the (smooth) solid Sh and 
the (smooth) plane P^. As the distance between the two is h + Xe, we deduce the expected 
lower bound. 

With similar- arguments, one has the upper bound: 

^h<^ + 0{\lnm. 

Indeed, let us define this time 

pO := {z = 0}, F^ := {x, x^Sh, z> 0}. 

Let be the Stokes flow in the smooth domain Fjj^. As the distance between Sh and P^ is h, 
the drag force satisfies 



/ |V<p = ^ + 0(|ln(/i)|) 



Now, can be extended by zero below P^ and defines in this way an element of Ah- In 



paiticulai", Fh < fp^ |Vt 
Hence, our methodology allows to derive quickly the inequalities 

6^ + 0{\ln{h + Xe)\) < Fh <^ + 0(|ln(/i)|). 



h + Xe h 

Interestingly, these bounds ai^e satisfied for any regime of parameters e and h. In particular, it provides 
the right asymptotic when e and h are of the same order. 
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Nevertheless, when e <C /i, it is fair to notice that a multiscale analysis gives a much refined 
description of the drag. For the sake of completeness, we briefly present it here. It relies on an 
asymptotic expansion of the Stokes flow = n| with respect to e. This expansion has already been 
described in close contexts, for instance in [jlSi . and we only recall its main elements. To keep track 
of the e dependency, we write P'^ instead of P, instead of Fh. We denote again P° and F° their 
smooth counterparts. 

The basic idea is to build an approximate solution app(^) (I2.1I )- (I2.2I )- (I2.3I ). in the form of an 
expansion in powers of e: 



+ ... + e^Uiy(x) + C/;f^(x,2/,x/e) ) (4.1) 



Each term of this expansion has two parts: 

• A regular part u]^ = ti5j(x) which models the macroscopic variations of the solution. 

• A boundary layer correction C/^ = C/^(x,y,X), which accounts for the fast variations of the 
solution near the oscillating boundary. Hence, it depends on the macroscopic variables x, y, but 
also on the microscopic variable X = x/e. It is defined for all 

x,y €R^, X={X, Y, Z) such that Z > 'y{X, Y). 

Moreover, C/^ is periodic in X, Y (due to the periodicity of the rough bondary 7 in X, Y) and 
satisfies 

limUi{x,y,X,Y,Z) =0, as Z ^ +00 

Back to the original variable x, this last condition corresponds to a boundary layer of typical 
size e near the rough wall P. 

Accordingly, the conesponding pressure field should read 

Ph,appi^) = Phi^) + Ph{x,y,^le) + e{pi{^) + Pl{x,y,i,/e)^ + ... 

We remind here the derivation of the 0(1) and 0{e) terms, which are enough for our purpose. 
First, if we inject the above expansions in (I2.11 i- (I2.2I) - (I2.3I) . and let Z — ;> +oo,we obtain by standai^d 
manipulations that 

-A4 + Vp^ = and div 4 = in F^, ul\as^ = e^, <|po = 0. 

Thus, we recover as expected that the leading term of the expansion is the Stokes flow without rough- 
ness. We then extend u\ and by zero below P^, so that they are defined over the whole Ff^. Such 
extensions trivially satisfy the Stokes equation for z < 0, as well as the no-slip condition at P^ . More- 
over, the velocity is continuous across the plane {z = 0}. But there is a 0(1) jump in the normal 
derivative. This explains the introduction of a boundary layer corrector with amplitude 0{e). Indeed, 
its gradient has amplitude 0(1), and allows to correct this artificial jump. 

Let us introduce the following notations: 



Vh{x,y,X) 
Phix,y,X) 
Phix,y,X) 



ui{x,y,0) + Ul{x,y,X), Z>j{X,Y), 
P^{x,y,X)-pl{x,y,0), Z > 0, 
P^{x,y,X), 0>Z>j{X,Y). 
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Note that, following the expansion (14.11 ). u*^ (x, y , 0) + eV^ (x, x/e) should be an approximation 
of the whole flow in the boundary layer. Plugging (14.11 ) into the equations, we derive fomially the 
following Stokes system: 

-AxVh + VxPh = 0, Z> 7(X, Y),Z^O, 
Vyr-Vh = 0, Z>0, 0, 
Vh = 0, Z = j{X,Y). 

together with the jump conditions 

Vh\z=o+ - VhU=o- = 0, {dzVh - Phez) \z=o+ - {dzVh - Phez) \z=q~ = -dzul{x, y, 0). 



Again, we stress that these jump conditions ensure the smoothness of the whole flow across the artifi- 
cial boundary {z = 0}. Note that by the divergence-free condition dzU^ ^{x, y, 0) = 0, so that only 
the horizontal components of dzu1{x, y, 0) are non-zero. Let us also point out that the variables x, y 
are just parameters in the system. In other words, one has 



Vh{x, y, X) = V(X)5,4(a;, y, 0), Ph{x, y, X) = V{X) ■ y, 0) 

for some 3-by-3 matrix function V and some 3d vector V which satisfy the (matricial) Stokes system 

-AxV + VxV = 0, Z>j{X,Y), 

Vx-V = 0, Z>jiX,Y), (4.2) 
V = 0, Z = j{X,Y), 

together with 

V|z=o+ - V|z=o- = 0, {dzV-r0ez)\z=o+-{dzV-r0ez)\z=o~ =-{111) 

This system of pde's, depending only on X, with periodic boundary conditions in the horizontal 
variable X, Y, has been extensively studied. We remind the following proposition, extracted from 
US: 




Proposition 1 The solution V of system (14.21) converges exponentially at infinity, that is 

\V{X,Y,Z) -V°^\ < Ce-^^ 
for some constant 3-by-3 matrix V°° and some 5 > 0. Moreover, V°° is of the form 

V°° = 

for some symmetric positive definite 2-by-2 matrix B. 

The non-zero block B is sometimes called the mobility tensor, see lfT6l . We stress that B is symmetric 
and definite positive, so diagonalizable in an orthonormal basis with positive eigenvalues. This fact 
will be used below. 

Back to Vh, we obtain that 

{Vh,x,Vh,y){x,y,X) Bdz {ul^x,uly) {x,y,0), Vh,z{x,y,X) 0, as Z +oo. 
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As the boundary layer coiTection Uf^ should decay at infinity, we obtain the boundary condition for 
the macroscopic coiTection u\ at P^. That is 

«,,<^) = ^a44,,4j, ul, = o, atp°. 

Together with the Stokes equations 

-Aul + Vpl = and div nj^ = in F^, 
and the boundary condition at the solid surface 

this determines u\, and ends the derivation of the 0{e) term of the expansion. The next order terms 
solve the same kind of equations, with inhomogeneous data coming from lower order profile. 

In a second step, one can show rigorously that the approximate solution is close to the exact 

solution n|. Indeed, introducing the differences 

•= '^h,app~'^h^ ^rid Q •= Ph,app~Ph 

leads to 

-Av + Vq = Rl, divv = rl, v\as^ = ifl- 

with remainder terms Rf^, rf^ and ipf^. For instance, the boundary data (pf^ is due to the boundary 
layer temis C/*(x,x/e), that do not vanish at dSh- We stress that the assumption e <^ h is cmcial 
for these remainders to be small. First, the boundary layer coiTcctions decay exponentially over a 
typical lengthscale e. To make it exponentially small at dSh, one needs e <^ h. Moreover, all other 
remainder terms are small with respect to e, but diverging with respect to h. Very roughly, they behave 
like 0{{e/h)^) where is the number of temis in the expansion (14. 11 1. The diverging powers of h 
come from taking derivatives of the u\, which are singular with respect to h. Again, the smallness 
condition e ^ /i is necessary. 

From there, as the remainder terms are small, one can through energy estimates deduce the small- 
ness of V, that is w In particular, for e small enough compared to h, the drag force on dSh 
reads 

- L - ^'•") • - ^ L. i'-^ - ^ -^'") • - ^ °<^'^ 

Moreover, it is easily seen that the fields uf^ := u^ + eu\, pf^ := p^ + ep\ satisfy the Stokes equation 
in Fj^, together with the boundary conditions 

UhlaSh = (^z, and ulJpo = 0, {uh,x^^h,y) \po = {uh,x,'^h,y) 1^^°- 

By the axisymmetry of F^, Kg = u^Re for any horizontal rotation Rg. As B is symmetric definite 
positive, this allows us to assume, up to a change of oithonormal basis, that B is diagonal with positive 
coefficients /S^, (3y. Now, there are two ways to interpret the effect of roughness. 

• On one hand, one can write the latter boundary condition as a slip condition of Navier type: 

|po = £l3dz{ul^^,ul^y)\po + o{e). 
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This is the so-called phenomenon of apparent slip, see lITSl . In the isotropic case (3 := = Py, 

one can use the bounds on the drag force derived in section |3] In the regime s^^h 1, this 
yields: 

^<J'h<T^ c,C>0. (4.3) 
n n 

• On the other hand, the drag force reads 



dn 



Using again the symmetry properties of u^^, we obtain that 

- /^^ {^-P>)-'^ - l^i\9Aj' + \dA,y\'), P ■■= ^4^- (4.4) 

But this last expression can be seen as the drag force created by a Stokes flow u^, between the 
solid Sji and a shifted wall Pg := z = —e /3. Indeed, following |[T9l . we get 



f2 



Using that 

ul{x, y, 0) = y, -e/3) + ef3dzU^{x, y, -e(3) + o{e) = e/Sd^ulix, y, 0) + o(e) 



we recover the same expression as in (I4.4l i. This interpretation of the roughness effect as a shift 
of the smooth wall yields 

-F. - (4.5) 

h + ep 

Note that bounds (14.31) and (14.51 ) are coherent, since we are here within the asymptotics e ^ /i ^ 1. Of 
course, these bounds would lead to very different behaviours for smaller h, as (14.31 ) forbids collision 
whereas (14.51) allows it. 



A Asymptotics of I and J 

In this appendix, we detail the computation of X(/3) and J'{(3, h) defined in (13.71) depending on the 
values of j3. 
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Case /3 < < 1 . When /3 is small, we expand with respect to /3 : 



„3 „3 „4+« / „5+2a 



This yields: 

/O (1 + ^)3 (1 + ^)4 



where routine calculations yield: 



Replacing in X(/3), we obtain, 

Case /? >> 1. When /3 is large, we split X(/3) = X0(/3) + X°°(/3) where: 



(l + ^ + /3si+")3' A (i + ^ + /3si+-)3' 

1 2 

To compute X°°{I3), we set s = j3'^~°'S and expand the integrand with respect to l/(/3i-"). This 
yields: 

Consequently, we distinguish three cases : 



• fora > 1/3, X°^(/3) = O ( ^ 



for Q = 1/3, X°°(/?) = + O ( ^ 



2/33 ' ^ ypz^ 
/i„ ^ / 1 \ s^ds 



fora<l/3, X-(/3) = -^ + 0(^-l) , /^a ^ ,.2 



;3T^ V/3V' Jo (f + §i+")3 

InX''(/3), we set u = /3s^+° and expand the integrand w.r.t. l//?i+^. This yields: 



(l + a)/3^ JO (1 + -")^ V/? 

with : 

. for a < 1/3, X0(/3) = O 
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.for a = 1/3, ^V) = ^^ + 0(^) 



.fora>l/3, X0(/3) = ^ + of-i), Tj^^ 

We obtain ( I3.9I ) comparing the values of and X°°(/3) in the three cases a<l/3, a = l/3 and 
Q > 1/3. 

It remains to handle the remainder term J'{f3,h). As previously, we split it into J'^{l3,h) + 
J°°{l3,h) where: 

/■I „7j„ rro/Vh „7j„ 



/o (1 + f + /3si+")4 7i (1 + ^ + /3si+")4 

1 

We set s = /3 1-°= s in the last integral. This yields: 

r.ro/(v^/3i/(i-)) ^7^g r 0(|ln(/i)|) if/3<l 



+ I 0(|ln(/3) + i^ln/i|) if/3>l. 
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